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We have obtained accurate ab initio 4 E + quartet potentials for the diatomic metastable triplet 
helium + alkali-metal (Li, Na, K, R.b) systems, using all-electron restricted open-shell coupled 
cluster singles and doubles with noniterative triples corrections [CCSD(T)] calculations and accurate 
calculations of the long-range Cq coefficients. These potentials provide accurate ab initio quartet 


scattering lengths, which for these many-electron 
masses and shallow potentials that results in a 
relevant for ultracold metastable triplet helium + 

PACS numbers: 31.15.A-, 34.20. Cf, 34.50. Cx, 67.85.-d 

I. INTRODUCTION 

Interactions and collisions involving helium in the 
metastable triplet 2 3 Si state (denoted as He*) have been 
regarded for many years as one of the most fascinating 
in atomic and molecular physics. The enormous amount 
of internal energy (19.8 eV) allows for Penning ioniza¬ 
tion, which has been extensively exploited in crossed- 
beam studies [I]. More recently, molecular-beam tech¬ 
niques involving He* have achieved sufficient resolution 
in kinetic energy M to observe for shape resonances in 
single partial waves for sub-Kelvin collisions, which are 
a sensitive probe of the interaction potential. 

Elastic and inelastic collisions at sub-milliKelvin ener¬ 
gies are relevant for ultracold trapped He* gases [5]. Pen¬ 
ning ionization greatly limits the lifetime of the trapped 
gas. However, for a spin-polarized gas of He* Penning 
ionization is strongly suppressed due to spin conservation 
[s]. Bose-Einstein condensates of 4 He* have been real¬ 
ized , benefiting from a sufficiently large scattering 

length of 142 ao [HI that allows for efficient evaporative 
cooling. Degenerate Fermi gases of 3 He* are obtained 
by sympathetic cooling with 4 He* [k|, which is efficient 
due to a very large interspecies scattering length of 496 a o 

Ell- 

One of the unique features of He* is its simple elec¬ 
tronic structure. Thus, it is being explored extensively 
in atomic physics to confront the state-of-the-art, ultra- 
precise spectroscopy and the most advanced fundamental 
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systems is possible, because of the small reduced 
small amount of bound states. Our results are 
alkali-metal mixture experiments. 


theories. For example, precise measurement of transition 
frequencies and the lifetimes are used to test the quantum 
electrodynamic (QED) calculations Il3 - (l6| . Similarly, 
the two-body interaction potential for spin-polarized He* 
atoms is so far the only system in ultracold physics for 
which it is possible to calculate the aforementioned scat¬ 
tering lengths |121 with an accuracy surpassing the exper¬ 
imental value [17| . This contrasts with other systems, 
such as e. g. the interaction potential for alkali-metal 
atoms, for which predicting the scattering length with 
such accuracy is impossible, and without experimental 
data, it is essentially unknown. 

Recently we have challenged this situation for the in¬ 
teraction between He* and Rb, and demonstrated a very 
good agreement between theoretical predictions based 
on all-electron restricted open-shell coupled cluster sin¬ 
gles and doubles with noniterative triples corrections 
[CCSD(T)] calculations and the scattering length derived 
from thermalization measurements of an ultracold mix¬ 
ture of 4 He* and 87 Rb [HI- The fact that even for so 
many electrons it is possible for ab initio calculations to 
quantitatively predict scattering lengths results not only 
from the simple electronic structure of He*, but mostly 
comes from the small reduced mass and shallow interac¬ 
tion potential that give a small number of bound states, 
making the scattering length less sensitive to uncertainty 
in the calculated potential (see Ref. 0 for the opposite 
case). 

Motivated by our previous work on He*+Rb [l8|, we 
consider the relevant interaction potentials for He* and 
the other alkali-metal atoms. Similarly as in the case of 
homonuclear He* collisions, Penning ionization: 


He* + A -a- He + A+ + e 


(1) 
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is suppressed by spin-polarizing He* and alkali-metal 
atom A both in the spin-stretched states, such that the 
total spin and its projection is maximum [20|. While 
the s character of the valence electron of He* and alkali- 
metal atoms are alike, we expect the amount of suppres¬ 
sion to be similar to that in He*+He* (in contrast to the 
other metastable noble gas systems |5|). Here we have 
assumed that the product ion A + has zero spin, which 
is true for the ground state. However, if excited, non¬ 
zero spin, (A + )* states are energetically available, Pen¬ 
ning ionization is not spin forbidden, even for the dou¬ 
bly spin-stretched state combination. Among the alkali- 
metal atoms the (A+)* states cannot be reached, except 
for Cs, for which the excitation energy from the neutral to 
the first excited ionic state is 17.2 eV (2lj |. Therefore we 
will discard Cs in this work as a stable ultracold mixture 
with He* seems experimentally not feasible and computa¬ 
tionally the presence of the energetically available excited 
ionic channel will complicate the calculations. Note that 
the above mentioned criteria to suppress Penning ioniza¬ 
tion exclude most other atomic species as well (except 
He*+H). 

The collision properties for He* + alkali-metal atom 
are determined by a 2 E+ doublet and a 4 E+ quartet 
potential, however, for the doubly-spin stretched state 
combination scattering only occurs in quartet potential. 
Therefore for realizing a stable mixture the properties of 
the quartet potential are the most relevant, in particular 
the quartet scattering length. This scattering length de¬ 
termines interspecies thermalization rates and therefore 
whether sympathetic cooling is efficient or not. Also, 
for quantum degenerate mixtures it determines, together 
with the intraspecies scattering lengths, whether the mix¬ 
ture is miscible or immiscible in case of Bose-Bose mix¬ 
tures [HI, [ 23 } , or whether a Fermi core or shell is formed 
in case of Bose-Fermi mixtures [24[ . Similar to the triplet 
potential in alkali-metal + alkali-metal interactions, the 
quartet potential is quite shallow, which in combination 
with the small reduced mass leads to a small number of 
bound states in the range of 11 to 15. Note that pre¬ 
vious experimental work on these kind of collisions sys¬ 
tems was based on measurin g el ectron emission spectra 
in crossed beam experiments |25l - l27j] , which inherently is 
only sensitive to the doublet potential. Therefore ultra¬ 
cold mixture experiments are the first to explore these 
quartet potentials. 

In this paper we present ab initio calculations of the 
quartet potentials of He* + alkali-metal (Li, Na, K and 
Rb) systems, where the results for He*+Rb are taken 
from Ref. M- In Sec. EH we describe the methodology 
of the calculations, while in Sec. m we present the ob¬ 
tained potentials and give a detailed discussion on the 
accuracy of those potentials. In Sec. IIVI we provide the 
corresponding scattering lengths for all the isotopologues 
and discuss possible implications for future experiments. 
Finally, in Sec. [V] we conclude and give an outlook. 

Throughout this paper we use Bohr radius, ao = 
5.2917721 x 10 -11 m, as a length unit and cm -1 = 


1.9864455 x 10 23 J as an energy unit. 


II. THEORY 


A. General considerations 


The s-wave scattering length a is obtained by solving 
the ID radial Schrodinger equation with zero angular mo¬ 
mentum and vanishing kinetic energy E: 

f(r) + |[£-y(r)]^(r) = 0, (2) 


where /r is the reduced mass, r the internuclear dis¬ 
tance and V (r) is the interaction potential. Current 
quantum chemistry ab initio methods are able to de¬ 
termine the short-range part of V(r) with an accuracy 
on the order of few percent, which translates into few 
cm -1 for dispersion-bound systems, like spin-stretched 
alkali-metal dimers or He*+alkali-metal atom systems. 
The long-range part of V (r) for alkali-dimers and sys¬ 
tems with helium atom can be described very accurately 
through the van der Waals expansion: the correspond¬ 
ing coefficients for such systems can be obtained with 
sub-percent accuracy. 

With an appropriate analytical form of V[r) it is 
possible to explore separately the influence of short- 
and long-range potential modifications on the scattering 
length. A good choice of such a function is the so-called 
Morse/Long-Range (MLR) potential [28}, which has the 
form: 


V(r) = D e 1 - 


MLR.(r) 

ULn(,r e ) 


exp [-</>(r)y p (r)} ) - D e , (3) 


where 

mlr(l) 

Vk(r ) 
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( 5 ) 

( 6 ) 


The free parameters of the MLR potential, determined 
by fitting, are the (f>j (j = 0,...,4) coefficients, while 
the potential well depth D e , equilibrium distance r e and 
4>oc = log[2D e /'ULR.('Ce)] are directly obtained from the ab 
initio calculations of the short-range potential. For the 
long-range part of the potential, V(r) r —>° —ulr(t) = 
—Cer~ 6 —Csr~ s —C ior~ 10 . Note that the statistical error 
introduced by the analytical fit is much smaller than the 
systematic uncertainty in the ab initio calculations. 


B. Short-range potential 

For the short-range potential we have used the cou¬ 
pled cluster method [29} in which the correlated elec- 
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tronic wavefunction is represented by the exponential op¬ 
erator exp(T) acting on the Slater determinant, where 
T = T\ + T 2 + ■ ■ . is the so-called cluster operator 
which includes single-, double- and higher-order excita¬ 
tions. A gold standard for weakly bound systems has be¬ 
come the approximate coupled cluster method denoted as 
CCSD(T), in which one fully includes single- and double¬ 
excitations and treats triple-excitations approximately. 
In this paper we use the open-shell version of CCSD(T) 
introduced by Knowles et al. [30], implemented in the 
MOLPRO program (31]. 

The accuracy of CCSD(T) for predicting binding en¬ 
ergies for weakly bound systems is typically on the or¬ 
der of 2-3% percent [32:]. Here we confirm this accu¬ 
racy by calculating the D e parameter using an even more 
sophisticated method, namely coupled cluster with full 
triple excitations, CCSDT. In addition, we have incorpo¬ 
rated scalar relativistic effects using the Douglas-Kroll- 
Hess approximation up to the fifth order in external po¬ 
tential (DKH5) |33]. All electrons have been correlated 
in these calculations. Interaction energies were calcu¬ 
lated using the counterpoise correction scheme of Boys 
and Bernardi [34j . in which the total monomer energies 
calculated in dimer basis sets are substracted from total 
energy of dimer. 

The accuracy of employed quantum chemistry method¬ 
ology is also determined by the choice of appropriate 
gaussian basis set describing the orbitals. Due to com¬ 
putational limitations, the basis set commonly used in 
quantum chemical calculations are truncated. However, 
the basis sets that are used in this paper are tailored 
to systematically improve the correlation energy and al¬ 
low approximate extrapolation to the complete basis set 
limit (CBS). By increasing the maximum angular mo¬ 
mentum in the basis sets the errors in correlation energy 
are expected to decrease. Hence, in order to converge 
the quantum chemical calculations one should obtain the 
results for a series of basis sets with increasing maximum 
angular momentum. 

While for most atoms corresponding to the first- and 
second rows in periodic table a variety of various fam¬ 
ilies of well-optimized, systematically convergent basis 
sets are available, for heavy alkali-metal atoms (K, Rb, 
Cs) the choice of basis sets is limited. We have used a 
sequence of core-valence correlation consistent basis sets 
developed by Prascher et al. [35] for Li and Na, which 
we will denote as TZ, QZ and 5Z. For K we have used 
an uncontracted atomic natural orbital relativistic basis 
sets [13 (denoted as ANO-RCC) with the exponents of 
h functions taken from g of the same basis (we follow the 
usual convention to label the gaussian functions with ap¬ 
propriate orbital angular momentum). In each case the 
basis set have been augmented by a set of two additional 
diffused functions per shell generated as an even tem¬ 
pered set according to prescription implemented in the 
MOLPRO program. The angular momentum structure of 
largest basis sets for Li, Na and K is 20slApl0d8f6glh 1 
22s20pl0d8f6g2h, 23sl8p7d4g2h, respectively. 


A proper choice of helium basis set is also essential: 
the optimal basis sets for the metastable triplet state 
have entirely different character compared to basis sets 
for ground state helium 0. Hence, we have decided 
to build a new basis set that takes into account diffused 
character of He*. To this end we have optimized a new 
set of exponents according to a following procedure. The 
starting point was an uncontracted ANO-RCC basis set 
for ground-state helium, i. e. 9s4p3d2f. Exponents were 
re-optimized to minimize the total energy in the DKH5 
relativistic method (at the full configuration interaction 
level) and extended to 15s8p5d4/2g, which gives less than 
1 /iHartree convergence. Then the lowest exponents for 
each shell were again augmented and re-optimized for the 
total energy of helium dimer quintet state at the equilib¬ 
rium distance. The final helium basis set has an angular 
momentum structure of 17sl0p7c?5/4g. 

With such basis sets we have performed test calcu¬ 
lations for the spin-polarized homonuclear systems He?!, 
Li 2 , Na 2 and K 2 to assess their performance near the 
equilibrium distance. We have obtained very good agree¬ 
ment of D e parameters with experimental results: for 
He^ we have obtained 1042.3 cm -1 which is very close 
to estimated limit of complete basis set for CCSD(T) 
method (1042.9 cm -1 ) [l2| and about 5 cm -1 shal¬ 
lower than the theoretical potential obtained with a 
full configuration interaction method. For Li 2 , Na 2 
and K 2 systems CCSD(T) values calculated with ba¬ 
sis sets described above are respectively, 330.5 cm -1 , 
172.3 cm -1 and 247.5 cm -1 , which can be compared with 
experimental values of 333.69 cm -1 , 174.96 cm -1 and 
255.017 cm -1 [37l - l39i ] . Clearly CCSD(T) with current 
basis sets systematically underestimates the well depths, 
but the deviation from the benchmark values are small. 

C. Long-range potential 

The Cq van der Waals coefficients of a general A+B 
system can be obtained by integration of the dipolar dy¬ 
namic polarizabilities a over the imaginary frequencies 

0 : 

c q B = — J a A (ico)a B (iu})du}. (7) 

For He* we have used the polarizabilities that are ob¬ 
tained from explicitly correlated gaussian wavefunctions 
with an accuracy on the order of 0 . 1 % for the zero- 
frequency [l 8 |, 0]. In case of alkali-metal atoms we have 
used the dynamic polarizabilities at imaginary frequen¬ 
cies given by Derevianko et al. 0]. These polarizabili¬ 
ties give homonuclear Cq coefficients with an accuracy of 
0.14%, 0.25% and 0.4% for the Li+Li, Na+Na and K+K 
systems, respectively. Hence, the inaccuracy of Cq in 

heteronuclear systems determined as \{5\ + £ 3 ) 0 

(where <5 a and <5 b are unsigned errors of pertinent static 
polarizabilities of monomers A and B, respectively) is al¬ 
ways smaller than 0.25%. Therefore the uncertainty in 
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FIG. 1. (Color online) Results of the ab initio calculations 
of the potential energy curves of the quartet ( 4 E + ) He* + 
alkali-metal systems. 


the long-range part of the potential is much smaller than 
that of the short-range part obtained from CCSD(T), 
and one can treat the long-range part of the interaction 
as fixed by theory. Note that our C$ coefficients agree 
with the values obtained by Zhang et al. [43| to better 
than 0.1%, which suggests that their estimated uncer¬ 
tainty of 1-5% is too conservative. 

For the Cg and C io coefficients we have used the values 
calculated by Zhang et al. f43| for which the authors 
estimated an uncertainty of about 1-10%. However, our 
study of He*+Rb system has shown that the error bounds 
for the C$ coefficient given by Zhang et al. was also far 
too conservative [l8j . 


III. AB INITIO POTENTIALS 
A. Recommended ab initio potentials 

The results of our calculations for He*+(Li, Na, K, Rb) 
quartet potentials are shown in the Fig. [T] the values of 
the potential well depth D e , the equilibrium distance r e 
and the Ce coefficient are given in Table Q] A complete 
list of the parameter values of the MLR potential are 
given in Appendix [A] (Table IIVI) . The uncertainty of the 
ab initio calculations is predominantly translated in an 
uncertainty in D e , while the uncertainty in r e is smaller 
than 0.01 a o- 

The overall qualitative pattern of the potentials resem¬ 
bles the triplet states of homonuclear alkali-metal sys¬ 
tems: the one containing Li has the deepest well, K and 
Rb are slightly shallower than Li, and Na is anomaly 
shallower than the other potentials. This pattern can be 
explained by subtle interplay between attractive disper¬ 
sion forces and Pauli repulsion (exchange energy). The 
Pauli repulsion increases with system size as the atoms 
are systematically more diffused and have systematically 


TABLE I. Key parameters of potential energy curves of quar¬ 
tet ( 4 S + ) He* + alkali-metal systems, including equilibrium 
distance r e , the potential depth D e and the Ce long-range 
coefficient. The uncertainty of the CCSD(T) calculations is 
reflected in the error bars on D e (see text). 


system 

r e (a 0 ) 

D e (cm 4 ) 

Ce (cm 1 ao) 

He* + Li 

7.53 

575±t 

4.5782 x 10 8 

He* + Na 

8.57 

36ll 4 

4.7723 x 10 8 

He* + I< 

9.08 

4701? 

7.7314 x 10 8 

He* + Rb 

9.41 

4531® 

8.4673 x 10 8 


larger radius. Increase in exchange energy is reflected 
in monotonic increase of equilibrium distance of dimers. 
On the other hand, the dispersion interaction gives rise 
to the attraction of atoms in the system and its magni¬ 
tude correlates with the C§ coefficients. Li and Na have 
comparable dispersion interactions, but a much larger 
exchange repulsion in the case of Na results in a much 
smaller well depth. On the other hand, K and Rb ex¬ 
hibit much stronger dispersion interaction, thus their well 
depths are noticeably larger than that of Na. 


B. Accuracy of the ab initio potentials 

To provide tests of the potential accuracy in the min¬ 
imum region we have performed additional calculations 
for He*+(Li, Na) usin g a coupled-clusters approach with 
full triple excitations [44| (CCSDT) for basis sets with 
maximum angular momentum limited to / (for Li and 
Na) and d (for He*) functions. CCSDT calculations were 
performed using the CFOUR quantum chemistry code [45| . 
The results are shown in the Table El and are discussed 
below. Note that a detailed description of the accuracy 
of the He*+Rb potential is given in Ref. fl8j . 


1. Basis set convergence 

We have studied the basis set convergence for all sys¬ 
tems. To this end we studied the dependence of the D e 
parameter in family of basis sets obtained by taking out 
from our actual basis sets two- and one- highest angular 
momentum functions. Such basis sets roughly correspond 
to triple- and quadruple-zeta quality (which we denote as 
TZ and QZ, respectively), while our best basis set is of 
5-zeta quality (5Z). 

We have found that even for the smallest basis sets 
of TZ quality the binding energies are very close to the 
values obtained from those obtained with 5Z basis sets: 
for the He*+(Li, Na) systems the CCSD(T) interaction 
energies calculated for the recommended r e ’s (Tab. HVl) 
are, respectively 569.2 cm -1 and 350.2 cm -1 , respec¬ 
tively, that is 6.2 -1 and 10.6 cm -1 below the recom¬ 
mended D e parameters. Such rapid convergence might 
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TABLE II. Interaction energies of quartet states of the 
He* + (Li, Na) systems calculated for the distance correspond¬ 
ing to r e of recommended potential (7.53 a o and 8.57 do, 
respectively) using CCSD(T) and CCSDT levels of coupled- 
clusters theory. We report also test calculations with various 
basis sets and levels of frozen-core approximation. See Sec. 
Ill Bl for notations regarding the basis sets. 


basis set 

level 

active electrons 

V(r e ) (cm J ) 

He*±Li 

TZ 

CCSDT 

all 

572.6 

TZ 

CCSD(T) 

all 

569.2 

TZ 

CCSDT 

3 

578.3 

TZ 

CCSD(T) 

3 

576.6 

QZ 

CCSD(T) 

all 

574.5 

5Z 

CCSD(T) 

all 

575.4 

CBS 

CCSD(T) 

all 

576.1 

He*±Na 

TZ 

CCSDT 

7 

352.0 

TZ 

CCSD(T) 

7 

350.2 

QZ 

CCSD(T) 

all 

362.2 

5Z 

CCSD(T) 

all 

360.8 

CBS 

CCSD(T) 

all 

359.3 


result from the fact that in the high-spin He*+alkali- 
metal systems same-spin electronic pairs are the main 
contribution to the interaction energy and it is known 
that the correlation energy for such pairs saturates faster 
than for the opposite-spin pairs. While we have esti¬ 
mated the CCSD(T) complete basis set (CBS) limit for 
the He*+(Li, Na) interaction energies to be 576.1 cm -1 
and 359.3 cm -1 , respectively, we rather prefer to treat 
the difference between 5Z basis set and CBS interaction 
energy as (unsigned) uncertainty attributed to the basis 
set incompleteness and to take the results for 5Z basis set 
as recommended values with basis set errors of ±0.6 cm -1 
and ±1.3 cm -1 , respectively. 

Our basis set for He*±K is not correlation consistent, 
so instead of performing the CBS estimate we have sim¬ 
ply compared the result in extended ANO-RCC basis set 
(see Sec. Ill Bl) to the interaction energy obtained with 
basis set with removed h functions. The value of the lat¬ 
ter is 468.6 cm -1 , hence the uncertainty due to basis set 
incompleteness is about ±1.3 cm -1 . 


2. Post-CCSD(T) contributions to the interaction energy 

Using the CCSDT method in reduced basis sets we 
have explored the performance of ab initio methods be¬ 
yond the CCSD(T) model for He*±(Li, Na) systems us¬ 
ing the TZ-quality basis sets. For the systems containing 
non-polar species quadruply excited configurations give 
substantially smaller contribution to the binding ener¬ 
gies [13, hence we might treat CCSDT interaction en¬ 


ergy at equilibrium as a probe of post-CCSD(T) effects 
near the equilibrium distance. Moreover, for the He*±Li 
case CCSDT is exact if the Is electrons of Li are kept 
frozen, and is exact also for isolated Li and He* atoms. 
We have calculated the CCSDT interaction energy for 
all electrons active and for the case when Is orbital is 
frozen. The results are given in the Table [Tl] In a basis 
set of TZ-quality the CCSDT interaction energy is deeper 
by merely 3.4 cm -1 . The TZ basis set is remarkably close 
to the complete basis set convergence limit (CBS), and 
we can safely assume that the difference between CCSDT 
and CCSD(T) is also nearly converged. It is interesting 
to notice that CCSDT for 3-electron calculations (i. e. 
with frozen Is shell of Li) gives actually an even smaller 
difference between CCSDT and CCSD(T) interaction en¬ 
ergies (1.7 cm -1 ), which shows that the core-relaxation 
effects in this case are comparable to the contributions 
beyond the CCSD(T) model. 

For He*±Na we were able to calculate the CCSDT 
interaction energies for 7 active electrons. It turns out 
that the difference between interaction energies CCSDT 
and CCSD(T) is only 1.8 cm -1 in TZ-quality basis set. 
We might assume that with converged basis set and all 
electrons correlated the error might be at most twice as 
large. 

For He*±K we were not able to converge the CCSDT 
calculations. Hence, to estimate bounds on D e we have 
compared how the analytical long-range potential given 
by van der Waals series compares to the interaction en¬ 
ergies from the actual calculations. It turns out that Cq 
coefficient extracted from the ab initio CCSD(T) inter¬ 
action energies (fitted from 20 ao to 35 ao) is about 1.5% 
smaller compared to the value obtained from perturba¬ 
tion theory (i. e. Eq. 0. If this value is treated as an 
estimate of ab initio potential uncertainty it translates 
to about ±7 cm -1 . 


3. Bounds on D e parameters 

As mentioned in Sec. Ill Bl our methodology predicts 
that the well depths of the homonuclear dimers He?;, LU, 
Na 2 , K 2 are systematically shallower compared to the 
experimental values by 0.53%, 0.92%, 1.52% and 2.94%, 
respectively. Because of the simple structure of the He* 
atom (for isolated He* atom CCSD(T) method is exact) 
we expect that within the He*±alkali-metal atom sys¬ 
tems the errors should be even smaller. By taking an 
average of the appropriate percentage uncertainties for 
homonuclear alkali dimers, i. e. (5ab = (<5a ± <5b)/2 we 
might expect that our ab initio potentials are deeper by 
about 0.7%, 1% and 1.7%, which translates to 4.0 cm -1 , 
3.7 cm -1 , 8.1 cm -1 for the He*±(Li, Na, K) potentials, 
respectively. These uncertainties are consistent with our 
estimate of post-CCSD(T) interaction energies, which in 
each case predicts small and systematically positive con¬ 
tributions beyond CCSD(T). We can conservatively as¬ 
sume that the real potentials have a well depth parameter 
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TABLE III. Scattering lengths for all He* + alkali-metal iso- 
topologues, showing the scattering length a corresponding to 
the recommended D e , and the bounds [a_;a + ] corresponding 
to the bounds on D e . Also the number of bound states N 
is given. Note that for the alkali-metal atoms the even iso¬ 
topes are fermions and the odd ones are bosons, while 4 He is 
a boson and 3 He is a fermion. 


system 

isotopes 

a 

[o; a+] 

N 

He*+Li 

3 + 6 

+26 

[+23; +26] 

11 


3 + 7 

-17 

[-27;-15] 

11 


4 + 6 

+22 

[+19; +23] 

12 


4 + 7 

-193 

[-607; -161] 

12 

He*+Na 

3 + 23 

+58 

[+52; +59] 

11 


4 + 23 

+7 

[-2; +9] 

12 

He*+K 

3 + 39 

+51 

[+42; +52] 

13 


3 + 40 

+49 

[+41;+51] 

13 


3 + 41 

+48 

[+39; +49] 

13 


4 + 39 

+97 

[+74; +101] 

15 


4 + 40 

+91 

[+70; +94] 

15 


4 + 41 

+86 

[+67; +89] 

15 

IIe*+RlW 

3 + 85 

+5 

[-17;+7] 

13 


3 + 87 

+3 

[-19;+5] 

13 


4 + 85 

+16 

[-4;+18] 

15 


4 + 87 

+15 

[-6;+17] 

15 


a The scattering length values are shifted by about 1 ao from our 
earlier reported values in Ref. list , which suffered from a small 
computational error. 


D e that is: i) not smaller than D e from CCSD(T) mi¬ 
nus basis set uncertainty; ii) not larger than our D e plus 
basis set uncertainty, plus the post-CCSD(T) correction. 
Hence, for He*+(Li, Na, K) systems the D e parameters 
for quartet potentials have uncertainties of -1/+4 cm -1 , 
-1/+4 cm -1 and 1/+9 cm -1 , respectively. 

Finally, we note that nonadiabatic effects can be ne¬ 
glected here despite the relatively small reduced masses. 
For the He*+Li system in the van der Waals mini¬ 
mum the diagonal Born-Oppenheimer correction calcu¬ 
lated with the Hartree-Fock electronic wavefunctions (us¬ 
ing the CFOUR program S3) is about 0.02 cm 1 . For the 
He* + (Na, K) systems these errors will be even smaller. 


IV. SCATTERING LENGTHS 

With the recommended MLR potentials the quartet 
scattering length a for all isotopologues can be calculated, 
for which we have used the ID renormalized Numerov 
propagator j4|| for a kinetic energy of 10 nK. Within 
the Born-Oppenheimer approximation it is only the dif¬ 
ferent reduced masses that give rise to different scatter¬ 
ing lengths for the isotopologues within each system and 
we have verified that non-adiabatic terms are negligible. 
We can conveniently explore the bounds on the ab ini- 



FIG. 2. (Color online) Scattering length as function of re¬ 
duced mass for the He*+Li potential, indicating the reduced 
masses for the four isotopologues (red dashed lines). The gray 
area gives the bounds on the scattering length related to the 
uncertainty in the potential. 


tio scattering lengths solely by scaling the D e parameter 
within its estimated bounds. Such scaling does not affect 
the long-range part of the MLR potential, which is kept 
fixed. The results are presented in Table m showing 
a and the bounds [a_;a+] related to the bounds on the 
calculated potentials, where the lowest possible value of 
D e corresponds to a+ and the highest one to a_. Also 
the number of bound states N is indicated, which differs 
between the two He isotopes by one or two units. 

The sensitivity of the scattering length a for the poten¬ 
tial well depth D e is very non-linear, and a diverges at 
each value of D e at which the potential supports a new 
bound state. In most cases here a is far away from a pole 
and the bounds on a are quite tight. A noticeable excep¬ 
tion is 4 He*+ 7 Li, which lies very close to a pole, leading 
to a broader range of possible scattering length values. 
This is also illustrated in Fig. [2j showing a as function 
of the reduced mass fi for the He*+Li potential. The ap¬ 
pearance of poles in a when varying D e or /r is similar, as 
/i and V ( r ) appear only as a product in the Schrodinger 
equation (see Eq. H). For He*+Li the four isotopologues 
have quite different reduced masses such that the corre¬ 
sponding scattering lengths can be very different. The 
similarity between the scattering lengths of 3 He*+ 6 Li 
and 4 He*+ 6 Li is purely accidental. For the heavier alkali- 
metal atoms (K and Rb) the scattering lengths for the 
same helium isotope are nearly the same, as the reduced 
mass hardly changes for the different alkali isotopes. 

The 16 isotopologues of He*+(Li, Na, K, Rb) contain 
6 Bose-Bose, 8 Bose-Fermi and 2 Fermi-Fermi mixtures. 
From the calculated scattering lengths we find that the 
4 He*+( 23 Na, 39 K, 41 K, 87 Rb) Bose-Bose mixtures are 
miscible, while the Bose-Fermi mixtures provide both 
Fermi core and Fermi shell situations. The 3 He*+ 40 K 
Fermi-Fermi mixture has already been proposed as its 
mass ratio is very close to a narrow interval where a 
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purely four-body Efimov effect is predicted [13 ■ It is 
interesting to note that both 3 He* and 40 K have an in¬ 
verted hyperfine structure, such that the high-field seek¬ 
ing doubly spin-stretched state combination is the lowest 
channel within the 3 He*+ 40 K manifold. 

From the calculated scattering lengths one can also 
find whether alkali-metal atoms can be used to sym¬ 
pathetically cool He*, as an alternative to evaporative 
cooling of 4 He* and sympathetic cooling of 3 He* by 
4 He*. Although both these schemes are successfully ap¬ 
plied, sympathetically cooling with another species would 
put a less stringent requirements on the initial num¬ 
ber of laser-cooled 4 He* atoms. Here one has to take 
into account that the thermalization rate scales with 
£ = 4mHe*mA/(mHe* +wa) 2 , and a 2 in the zero temper¬ 
ature limit (see for instance Ref. 0 )- Na might be the 
most suitable candidate for sympathetic cooling of 3 He*, 
while 'Li is the best coolant for 4 He* (although only until 
quantum degeneracy is reached because of immiscibility). 


V. CONCLUSIONS AND OUTLOOK 

We have obtained accurate ab initio 4 E + quartet po¬ 
tentials for He*+(Li, Na, K, Rb), using CCSD(T) calcu¬ 
lations and accurate calculations of the Cq coefficients, 
and have calculated the corresponding scattering lengths 
for all the isotopologues. An accurate prediction of scat¬ 
tering lengths for these many-electron systems is possi¬ 
ble, in contrast to nearly all other types of ultracold mix¬ 
tures, because of the small reduced masses and shallow 
potentials that results in a small amount of bound states 
(N = 11 — 15), and therefore a reduced the sensitivity of 
the scattering length to the properties of the potentials. 

So far, we have only considered the quartet potential, 
which is the only relevant potential for the doubly spin- 
stretched state combinations, for which Penning ioniza¬ 
tion is suppressed. Feshbach resonances that allow to 
tune the scattering length are in principle possible for 


He* + alkali-metal systems, as both atoms have elec¬ 
tron spin and at least one has nuclear spin, however they 
would require spin-state combinations in which at least 
one of the atom is not in the spin-stretched state, and 
scattering has both doublet and quartet character. Un¬ 
fortunately, accurate ab initio calculations of the dou¬ 
blet potentials are far more challenging. First of all, the 
doublet potentials need multiconfigurational treatment 
and these methods are at present far less accurate than 
CCSD(T). Secondly, for the doublet states coupling with 
continuum states of ionized channels is possible, which 
might complicate the calculations. In addition, the much 
larger well depths [26| , and therefore much larger amount 
of bound states, results in a much stronger sensitivity of 
the doublet scattering length to the underlying potential. 
Therefore, experimental input, such as positions of Fesh¬ 
bach resonances, will be needed to determine the doublet 
scattering properties. 
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Appendix A: Parameters of MLR potential 

In Table ITVl we give the parameter values for the MLR 
potentials, where D e and r e are obtained directly from 
the CCSD(T) calculations, C$ is calculated from the 
dynamical polarizabilities, Cg and Cio are taken from 
Ref. [0 , and the (f> parameters are obtained from fitting 
the MLR potential to the CCSD(T) data. The values of 
p and q are chosen to obtain the best fit. 


[1] P. E. Siska, Rev. Mod. Phys. 65, 337 (1993). 

[2] A. B. Henson, S. Gersten, Y. Shagam, J. Narevicius, and 
E. Narevicius, Science 338, 234 (2012). 

[3] E. Lavert-Ofir, Y. Shagam, B. Henson, Alon, S. Gersten, 
J. Klos, P. S. Zuchowski, J. Narevicius, and E. Narevi¬ 
cius, Nature Chemistry 6, 332 (2014). 

[4] J. Jankunas, K. Jachymski, M. Hapka, and A. Oster- 
walder, J. Chem. Phys. 142, 164305 (2015). 

[5] W. Vassen, C. Cohen-Tannoudji, M. Leduc, D. Boiron, 
C. Westbrook, A. Truscott, K. Baldwin, G. Birkl, P. Can- 
cio, and M. Trippenbach, Rev. Mod. Phys. 84, 175 
( 2012 ). 

[6] G. V. Shlyapnikov, J. T. M. Walraven, U. M. Rahmanov, 
and M. W. Reynolds, Phys. Rev. Lett. 73, 3247 (1994). 


[7] A. Robert, O. Sirjean, A. Browaeys, J. Poupard, 
S. Nowak, D. Boiron, C. I. Westbrook, and A. Aspect, 
Science 292, 461 (2001). 

[8] F. Pereira Dos Santos, J. Leonard, J. Wang, C. J. Bar- 
relet, F. Perales, E. Rasel, C. S. Unnikrishnan, M. Leduc, 
and C. Cohen-Tannoudji, Phys. Rev. Lett. 86, 3459 
( 2001 ). 

[9] A. S. Tychkov, T. Jeltes, J. M. McNamara, P. J. J. Tol, 
N. Herschbach, W. Hogervorst, and W. Vassen, Phys. 
Rev. A 73, 031603(R) (2006). 

[10] R. G. Dali and A. G. Truscott, Opt. Commun. 270, 255 
(2007). 

[11] M. Keller, M. Kotyrba, F. Leupold, M. Singh, M. Ebner, 
and A. Zeilinger, Phys. Rev. A 90, 063607 (2014). 



TABLE IV. Parameter values of the MLR potentials. The 
number of digits given exceeds their precision. The Cg and 
Cio coefficients are taken from Ref. 14311 . 


He*+Li 

D e 

575.40 cm -1 


4>o 

-2.5087 

r e 

7.5296 a 0 


<t> i 

0.32009 

Co 

4.5782x10 s cm -1 a® 


0 2 

-0.38608 

Cg 

2.9058X10 10 cm -1 a® 


0 3 

0.34763 

Cio 

2.8093xl0 12 cm -1 aj° 


0 4 

-0.78801 

P , Q 

4,4 


05 

-1.3009 

He*+Na 

D e 

360.80 cm 1 


00 

-1.8898 

r e 

8.5678 a 0 


01 

0.30840 

c 6 

4.7723x10 s cm -1 a® 


02 

-0.059153 

Cg 

3.4084X10 10 cm _1 a s 


03 

0.20557 

Cio 

3.4458x 10 12 cm _1 ao° 


04 

-0.86720 

P, Q 

5, 5 


05 

-0.47504 

He*+K 

D e 

469.83 cm" 1 


00 

-1.8233 

r e 

9.08326 a 0 


01 

0.33689 

C & 

7.7314x10 s cm -1 a® 


02 

-0.12801 

c 8 

6.7049X10 10 cm _1 a s 


03 

0.13803 

C\o 

7.6487X10 12 cm -1 aj° 


04 

-0.67352 

p , q 

5, 5 


05 

-0.22291 

He*+Rb 

D e 

452.71 cm -1 


00 

-1.8284 

r e 

9.4079 cio 


0i 

0.48678 

C 6 

8.4673x10 s cm -1 a® 


02 

-0.065081 

Cg 

8.0108X10 10 cm _1 a s 


03 

-0.30087 

Cio 

9.4242X 10 12 cm -1 ao° 


04 

-1.5195 

p, q 

5,4 





[12] M. Przybytek and B. Jeziorski, J. Chem. Phys. 123, 
134315 (2005). 

[13] J. M. McNamara, T. Jeltes, A. S. Tychkov, W. Hoger- 
vorst, and W. Vassen, Phys. Rev. Lett. 97, 080404 
(2006). 

[14] S. S. Hodgman, R. G. Dali, L. J. Byron, K. G. H. Bald¬ 
win, S. J. Buckman, and A. G. Truscott, Phys. Rev. 
Lett. 103, 053002 (2009). 

[15] R. van Rooij, J. S. Borbely, J. Simonet, M. D. Hooger- 
land, K. S. E. Eikema, R. A. Rozendaal, and W. Vassen, 
Science 333, 196 (2011). 

[16] R. P. M. J. W. Notermans and W. Vassen, Phys. Rev. 
Lett. 112, 253002 (2014). 

[17] S. Moal, M. Portier, J. Kim, J. Dugue, U. D. Rapol, 
M. Leduc, and C. Cohen-Tannoudji, Phys. Rev. Lett. 
96, 023203 (2006). 

[18] S. Knoop, P. S. Zuchowski, D. K§dziera, L. Mentel, 
M. Puchalski, H. P. Mishra, A. S. Flores, and W. Vassen, 
Phys. Rev. A 90, 022709 (2014). 


[19] M. Borkowski, P. S. Zuchowski, R. Ciurylo, P. S. Juli¬ 
enne, D. K§dziera, L. Mentel, P. Tecmer, F. Miinchow, 

C. Bruni, and A. Gorlitz, Phys. Rev. A 88, 052708 
(2013). 

[20] L. J. Byron, R. G. Dali, W. Rugway, and A. G. Truscott, 

New. J. Phys. 12, 013004 (2010). 

[21] J. E. Sansonetti, J. Phys. Chem. Ref. Data 38, 761 
(2009). 

[22] B. D. Esry, C. H. Greene, J. P. Burke, and J. L. Bohn, 

Phys. Rev. Lett. 78, 3594 (1997). 

[23] C. K. Law, H. Pu, N. P. Bigelow, and J. H. Eberly, Phys. 

Rev. Lett. 79, 3105 (1997). 

[24] K. Mplmer, Phys. Rev. Lett. 80, 1804 (1998). 

[25] L. G. Gray, R. S. Keiffer, J. M. Ratliff, F. B. Dunning, 
and G. K. Walters, Phys. Rev. A 32, 1348 (1985). 

[26] M.-W. Ruf, A. J. Yencha, and H. Hotop, Z. Phys. D 5, 

9 (1987). 

[27] A. Merz, M. Muller, M.-W. Ruf, and H. Hotop, Chem. 

Phys. 145, 219 (1990). 

[28] R. J. LeRoy, Y. Huang, and C. Jary, J. Chem. Phys. 

125, 164310 (2006). 

[29] J. Cizek, J. Chem. Phys. 45, 4256 (1966) 

[30] P. J. Knowles, C. Hampel, and H. J. Werner, J. Chem. 

Phys. 99, 5219 (1993). 

[31] H.-J. Werner, P. J. Knowles, G. Knizia, F. R. 

Manby, M. Schiitz, et al, “MOLPRO, version 2012.1, 
a package of ab initio programs,” (2012), see 
http: //www. molpro.net. 

[32] D. G. Smith, P. Jankowski, M. Slawik, H. A. Witek, 
and K. Patkowski, J. Chem. Theory Comput. 10, 3140 
(2014). 

[33] M. Reiher and A. Wolf, J. Chem. Phys. 121, 2037 (2004). 

[34] S. F. Boys and F. Bernardi, Mol. Phys. 19, 553 (1970). 

[35] B. P. Prascher, D. E. Woon, K. A. Peterson, T. H. Dun¬ 
ning, Jr., and A. K. Wilson, Theor. Chem. Acc. 128, 69 
( 2011 ). 

[36] B. O. Roos, R. Lindh, P.-A. Malmqvist, V. Veryazov, 
and P.-O. Widmark, J. Phys. Chem. A , 11431 (2008). 

[37] C. Linton, F. Martin, A. Ross, I. Russier, P. Crozet, 

A. Yiannopoulou, L. Li, and A. Lyyra, J. Mol. Spec- 
trosc. 196, 20 (1999). 

[38] T.-S. Ho, H. Rabitz, and G. Scoles, J. Chem. Phys. 112, 

6218 (2000). 

[39] A. Pashov, P. Popov, H. Knockel, and E. Tiemann, Eur. 

Phys. J. D. 46, 241 (2008). 

[40] P. Knowles and W. Meath, Mol. Phys. 60, 1143 (1987). 

[41] M. Puchalski (private communication). 

[42] A. Derevianko, S. G. Porsev, and J. F. Babb, At. Data 
Nucl. Data Tables 96 (2010). 

[43] J.-Y. Zhang, L.-Y. Tang, T.-Y. Shi, Z.-C. Yan, and 
U. Schwingenschlogl, Phys. Rev. A 86, 064701 (2012). 

[44] J. D. Watts and R. J. Bartlett, J. Chem. Phys. 93, 6104 
(1990). 

[45] J. F. Stanton, J. Gauss, M. E. 

Harding, P. G. Szalay, et al, 

“CFOUR, coupled-cluster techniques for computational chemistry, 
(2012), see http://www.molpro.net. 

[46] B. R. Johnson, J. Chem. Phys. 69, 4678 (1978). 

[47] Y. Castin, C. Mora, and L. Pricoupenko, Phys. Rev. 

Lett. 105, 223201 (2010). 

















